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Abstract 

Lattice  structures  exhibit  good  finite  word  length 
properties.  In  some  real-time  applications,  high  speed  is 
very  crucial  to  the  system  performance  along  with  finite 
word  length  effects.  Although  lattice  filters  can  be 
pipelined  by  cutset  localization  procedure,  no  sample 
rate  increase  can  be  achieved  by  that  method. 

In  this  paper,  a  pipelining  method  in  lattice  digital 
filters  is  introduced.  This  pipelining  method  is  based 
upon  constrained  IIR  digital  filter  design  method  by 
which  pipelined  direct-form  filters  are  designed.  These 
direct-form  filters  are  transformed  to  pipelined  lattice 
digital  filters.  It  is  shown  that  the  roundoff  error  and  the 
number  of  multiplytadd  operations  of  the  resulting 
pipelined  lattice  filters  are  smaller  than  those  of  the 
pipelined  lattice  filters  obtained  by  applying  look-ahead 
on  direct-form  nonpipelined  digital  filters. 


1.  Introduction 

Since  lattice  structures  exhibit  very  good  finite  word 
length  properties,  these  structures  are  preferable  for  fixed- 
point  finite  >vord-length  implementation.  Different  lattice 
digital  filter  structures  have  been  proposed  in  [1-4]. 

Achieving  high  speed  in  direct-form  and  and  lattice 
digital  filters  is  difficult  because  of  the  feedback  loop. 
Although  these  filters  can  be  pipelined  by  cut-set 
localization  procedure  proposed  in  [5],  it  should  be  noted 
that  no  sample  rate  increase  can  be  achieved  by  that 
method.  Cutset  transformation  can  increase  the  clock 
speed,  but  cannot  increase  the  sample  rate,  since  multiple 
clock  cycles  are  needed  to  process  one  sample. 

Implementation  of  high-speed,  pipelined,  stable, 
recursive  digital  filters  was  considered  impossible  before. 
Only  recently,  scattered  look-ahead  techniques  were 
proposed  to  convert  a  nonpipelined  filter  to  a  pipelinable 
form  [6].  The  scattered  look-ahead  technique  adds  poles  at 
equal  angular  and  radial  spacing  for  each  pole  of  the 
nonpipelined  filler,  and  these  poles  are  canceled  by  zeros. 
This  transformation  converts  the  denominator  of  the 
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transfer  function  to  a  pipelined  form.  This  technique,  for 
the  first  time,  made  it  feasible  to  implement  high  speed 
recursive  digital  filters  using  pipelining. 

One  drawback  of  the  scattered  look-ahead  technique  is 
the  introduction  of  the  canceling  zeros.  These  zeros 
increase  the  number  of  multiplication  operations  needed 
to  implement  the  digital  filter.  Therefore,  as  suggested  in 
[6],  it  is  of  interest  to  design  pipelinable  digital  filter 
transfer  functions  directly  from  the  filter  spectrum,  as 
opposed  to  designing  a  nonpipelined  filter  and  applying 
scattered  look-ahead  on  this  filter. 

To  this  end,  we  address  the  design  of  pipelined  direct- 
form  recursive  digital  filters  from  the  filter  spectrum.  By 
constraining  the  location  of  poles  at  equal  angular  and 
radial  spacing,  we  can  design  a  transfer  function  whose 
denominator  is  pipelinable.  Fortunately,  this  is  possible 
by  using  a  variation  of  the  filter  design  procedure  used  for 
decimation  filler  design  [7-10].  The  pipelinable  filter 
transfer  function  can  be  mapped  to  many  lattice  digital 
filter  structures.  Pipelining  the  orthogonal  structure  in 
[3]  leads  to  no  increase  in  speed.  Therefore,  the 
orthogonal  lattice  filter  structure  is  not  suitable  for 
pipelining.  The  minimum-noise  lattice  structure  in  [4] 
requires  more  number  of  muilliply-add  operations  than 
[1],  but  leads  to  less  roundoff  error.  In  ^is  paper,  we 
propose  the  use  of  pipelining  for  the  basic,  normalized, 
and  the  minimum-noise  lattice  suucture. 

This  paper  is  organized  as  follows.  The  Schur 
algorithm  [11]  which  forms  the  basis  for  design  of  lattice 
digital  fillers  is  reviewed  in  section  2.  Section  3  describes 
the  theory  of  lattice  filler  pipelining.  Section  4  presents 
the  constrained  filter  design  procedure  for  design  of 
pipelined  direct-form  recursive  digital  filters,  which  are 
then  reduced  to  lattice  structures.  Section  5  presents 
design  examples  for  several  lattice  filters. 
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2.  Schur  Algorithm 


One  of  the  important  properties  of  the  Schur 
algorithm  is  that  all  polynomials  expanded  from  any 
given  polynomial  by  the  Schur  algorithm  are  orthogonal 
to  each  other.  Schur  polynomial  is  a  polynomial  which 
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does  not  have  zeros  on  or  outside  the  unit  circle. 
Thereiv..e,  the  denominator  of  a  stable  IIR  filter  is  a 
Schur  polynomial.  Let  a  polynomial  with  real 
coefficients  be  defined  by 
N 

D(z)*]^diz‘. 

i=0 

Then,  initialize  ^th-order  Schur  polynomial  Ofi(z)  as 
<I>n(z)*D(z). 

N 

where  0^(2)  =  X 
isO 

From  4>|^z),  form  the  polynomial  4>i4_i(z)  as  follows: 


where  ^]^(z)  is  the  reverse  polynomial  of  <t>i4(z)  and  is 

defined  by  z^<I>n(z'’).  The  degree  of  «I>n_i(z)  is  1  less 
than  d>N(z). 

In  order  to  implement  the  algorithm,  each  of  the 
polynomials  d>N-i(z)  is  generated  from  <bN-i+t  (2).  for 
i=l,2, ...  in  the  same  manner  as  «1>n-i(2)  ‘s  generated 
from  d>N(z).  It  may  be  noted  that  the  coefficients  of 
increasing  powers  of  z  in  <I>n_i(z)  are  the  N  determinants 
of  2  X  2  submatrices  formed  by  the  first  column  and  each 
succeeding  column  in  the  following  matrix. 

J_  ^  ♦n-2  —  <1*1  1*0 

C 

<>0  1*1  ^2  •••  1*N-1  1*N  . 

where  C  =  V  1*N  *  1*0  • 

For  each  2x2  submatrix,  if  a  column  is  composed  of  all 
zero  elements,  the  determinant  of  the  corresponding 
submatrix  is  zero.  Using  this  property  of  determinant,  we 
can  observe  an  useful  fact  for  pipelining: 

If  OnCz)  is  a  Schur  polynomi^  of  order  N,  and  has  j- 

consecutive  zero  coefficients  between  each  two  nonzero 
coefficients  of  nearest  degree,  then  <I>N_i(z)  of  order  N-i 
also  has  j-consecutive  zero  coefficients  between  each  two 
nonzero  coefficients  of  nearest  degree  for  i  =  1  to  N-j-1 
and  has  only  one  nonzero  coefficient  for  i  =  N-j  to  N. 

3.  Pipelining  of  Recursive  Digital 
Lattice  Filters 

In  this  section,  we  study  pipelinability  of  basic. 


normalized  and  minimum-noise  lattice  filters. 
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Fig.  I  The  structure  of  Basic  laiiice  filter 

3.1  Pipelining  of  Basic  Lattice  Filter 

Fig.  1  shows  the  structure  of  the  basic  lattice  filter  which 
is  characterized  by  N  k-parameters  and  N+J  tap  gain 
parameters.  The  multiply-add  operations  in  the  feed¬ 
forward  section  can  be  pipelined  at  any  desired  level  by 
placing  latches  at  appropriate  feed-forward  cutset 
locations.  As  can  be  seen  from  Fig.  1,  the  feedback  loop 
in  a  module  is  removed  if  the  k-parameter  of  the  module 
is  zero.  Therefore,  by  removing  (M  -  /^-consecutive  k- 
parameters,  the  basic  lattice  filter  can  be  pipelined  by  Af- 
stages,  leading  to  M  times  increase  in  sp^  as  compared 
with  the  original  basic  lattice  filter. 

Let  the  A^th-ordcr  denominator  of  a  digital  filler  transfer 
function  be  denoted  by 

i=0 

where  m  -  N,  and  the  leading  coefficient,  d„  g,  is 

assumed  to  be  one.  Then,  k-parameters  are  recursively 
obtained  by  starting  firom  Dj»;(z)  using: 

D„(2)-k„..lDJz) 

2)D„.,(z)= - -2 - 

for  m=N.  N-1 . /  . 

As  can  be  seen  from  above  steps,  D„(z)’s  for  m  =  N, 
N-I.  ...  0  are  Schur  polynomials  with  a  different 
normalization  constant.  Therefore,  if  D^Cz)  has  (M  - 1)- 
consccutive  zero  coefficients  between  each  two  nonzero 
coefficients  of  nearest  degree,  then  D„(z)'s  for  m-N,N- 

1 . Af,  also  have  (M  -  /^consecutive  zero  coefficients 

between  each  two  nonzero  coefficients  of  nearest  degree 


and  D„(2)'s  for  m  =  M-1,  M-2 . 1.  have  only  one 

nonzero  coefficients.  Then,  since  k-parameters  are  the 

coefficients  of  z*“  of  each  D^(z)  for  m=N,  N-1 . 1, 

we  can  obtain  (M  •  inconsecutive  zero  k-parameters. 
which  leads  to  Af -stage  pipielining.  In  general,  an  M- 
stage  pipelined  filter  contains  Af -delays  in  every  loop. 
Since  the  number  of  multiply-add  operations  in  every 
loop  remains  constant,  the  Af-stage  pipelined  filter  can  be 
clocked  at  Af -times  faster  rate,  as  compared  with  the 
ncHipipelined  filter. 

3.2  Pipelining  of  Normalized  Lattice  Filter 

The  structure  of  normalized  lattice  filter  is  shown  in 
Fig.  2,  where  the  feedback  section  is  again  described  by 
k-paramelers  only.  The  k-parameters  are  calculated  by  the 
same  procedure  as  in  section  3.1. 

As  can  be  seen  from  Fig.  2,  there  are  four  multipliers 

within  a  module  nr.  k„,  -k„,  and  two  1  -  k^^  *5  jf 

is  zero,  then  1  -  k^^  becomes  one.  Therefore,  module 
m  needs  no  multipliers  and  the  feedback  loop  of  the 
module  is  removed,  which  means  the  normalize  lattice 
filter  can  be  pipelined. 
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Fig.  2  The  stiuclure  of  Noimalized  Uilice  filler 

3.3  Pipelining  of  Minimum-Noise  Lattice 
Filter 

The  structure  of  the  minimum-noise  lattice  filter  [4]  is 
shown  in  Fig.  3.  Feedback  section  is  described  by  k- 
parameters  only  and  the  k-parameters  are  calculated  using 
Schur  polynomials  as  follows: 


where  C>^(z)  is  mth-order  Schur  polynomial  obtained 


from  the  denominator.  If  k„  of  the  module  m  is  zero, 
then  the  feedback  loop  of  the  module  is  removed  since 
■^1  -  k„2  becomes  one.  For  k^  to  be  zero,  *^(0),  or 
the  constant  term  of  should  be  zero.  As  was 

discussed  in  section  2,  if  the  denominator  has  (M  - 1)- 
consecutive  zero  coefficients  between  each  two  nonzero 
coefficients  of  nearest  degree,  then  the  minimum-noise 
filter  can  be  pipelined. 


Fig.  3  The  sinicture  of  Minimum-noise  laiiice  filler 


4.  Design  of  Pipelined  HR  Lattice  Filter 
Using  Constrained  Filter  Design 

It  is  apparent  from  previous  sections  that  the  basic, 
normalized,  and  minimum-noise  lattice  filters  can  be 
pipelined  by  Af -stages  if  the  denominator  of  a  transfer 
function  has  fAf-7^<onsecutive  zero  coefficients  between 
each  two  nonzero  coefficients  of  nearest  degree.  The 
pipelinable  transfer  functions  can  be  obtained  by 
applying  the  scattered  look-ahead  method  to  the 
nonpipelined  filter  transfer  functions.  However,  by  using 
constrained  filter  design  method,  we  can  design  inherently 
pipelinable  transfer  functions  directly  from  the  filter 
spectrum  specifications. 

Our  procedure  is  a  variation  of  the  procedures  in  [7]  and 
[9].  The  design  method  in  [7]  first  expresses  the 
magnitude  and  group  delay  responses  of  the  filter  as 
functions  of  the  radii  and  angles  of  the  poles  and  zeros. 
Then  we  obtain  formulae  for  the  partial  derivatives  of 
the  magnitude  and  group  delay  with  respect  to  the  radius 
and  the  angle  of  a  pole  and  a  zero.  These  derivatives  are 
used  in  the  Fletcher-Powell  algorithm  to  minimize  the 
approximation  error  [8].  In  this  algorithm,  the  minimum 
value  is  obtained  by  only  using  the  first  derivative 
information.  To  obtain  an  inherently  A/-stage  pipelinable 
filter,  we  recompute  the  partial  derivatives  for  a 
denominator  in  powers  of  rather  than  2  [9]  and 
substitute  these  equations  into  the  program  in  [10]. 


Then,  the  resulted  denominator  is  in  terms  of  z*^. 

The  modified  filter  design  procedure  is  summarized  as 
follows: 

1)  For  the  given  filter  specifications  (pass-band,  stop- 
band.  pass-band  &  stop-band  ripples,  and  A/-stage 
pipelining),  we  start  with  1 -complex  pole  pair.  Then,  the 
order  of  the  denominator  is  2M  and  the  number  of  zeros 
are  restricted  so  that  the  order  of  the  numerator  is  less 
than  or  equal  to  2M.  We  usually  start  with  M  unit-circle 
zero  pairs. 

2)  If  the  filter  specifications  are  not  satisfied  by  1- 
complex  pole  pair,  we  increase  the  number  of  poles  and 
zeros  and/or  adjust  initial  positions  of  poles  and  zeros  and 
weighting  factors  for  the  pass-band  and  stop-band.  The 
niter  is  redesigned.  This  procedure  is  repeated  until  the 
filter  specifications  are  satisfied. 

3)  The  pipelined  direct-form  filter  transfer  function  is 
then  used  to  obtain  the  pipelined  lattice  filters  by  using 
the  sybthesis  procedures  of  section  3. 


5.  Design  Example  of  Pipelined  Lattice 
Filters 

Consider  the  lowpass  filter  with  I)pass-band:  0  0.2 
Jt.  ripple:  1  dB.  and  2)stop-band:  0.3  Jt  -  7C.  ripple:  15 
dB.  In  [12],  this  example  was  designed  by  using  the  4th 
order  Chebyshev  filter. 

N(z)  =  0.0018  +  0.0073z-‘  +  O.OOllz'^  +  0.0073z-3  + 
0.001 8z-4 

D(z)  =  1  -  3.0544z->  +  3.829Z-2  -  2.2925z-3  +  O.OSz-*, 
where  N(z)  and  D(z)  denote  numerator  and  denominator  of 
the  transfer  function,  respectively. 

Using  the  constrained  filter  design  procedure,  we 
obtained  the  uansfer  functions  for  different  number  of 
pipelining  stages.  These  filters  are  listed  below. 

1)  Nonpipelined 

N(z)  =  0.0307  +  0.0054z-»  +  0.0203z-2  +  0.0054z-3  + 
0.0307z-^ 

D(z)  =  1  -  2.5754z->  +  2.868z-2  -  1.5335r3  +  0.333z^. 

2)  2-stage  pipelined  (M  =  2) 

N(z)  =  0.0532  +  0.112Z-J  +  0.1159z-2  +  0.0941z-3  + 
0.1 159z-^  +  0.1  12z-5  +  0.0532z-« 

D(z)  =  1  -  1.02Z-2  +  1.065z-^  -  0.509z-«  +  0.128Z* 

3)  4-stage  pipelined  (M  =  4) 

N(z)  =  0.0754  +  0.132z-‘  +  0.2018z-2  +  0.241  IrU 
0.2815Z-'*  +  0.2411Z-5  +  0.201 8z-<5  +  0.1 32z-‘' + 
0.10754Z-*. 

D{z)  =  1  +  0.5626z-^  +  0.0791Z-*. 

4)  8-stage  pipelined  (M  =  8) 

N(z)  *  0.0974  +  0.0982z-‘  +  0.0814r2  +  0.0347z-3  + 
0.0047z  <  -  0.0275Z-5  +  0.029z^+  0.0761  z  '^  +  0.065z-* 
+  0.012Z-’  -  0.0434Z-10  -  0.058z-‘>  -  0.055z-‘2  - 
0.0034z>3  +  0.0082z-'<  +  0.0098z->5  +  0.0098z-‘6 
EKz) »  1  -  0.2013Z-*  +  0.101z*6 
Fig.  4  and  Fig.5  show  the  magnitude  response  and  the 


implementation  of  the  4-stage  pipelined  normalized 
lattice  filter  obtained  from  the  above  transfer  function. 
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Fig.  5  4'Siage  pipelined  Normalized  filler  by  Constrained  design 

The  feedback  loop  multiply-add  operations  can  be 
pipelined  by  redistributing  the  delays.  Roundoff  errors  of 
these  pipelined  normalized  lattice  filters  and  those 
obtain^  by  using  look-ahead  method  to  the  nonpipelined 
Chebyshev  filter  transfer  function  are  calculated  using  the 
method  introduced  in  [4],  and  are  compared  in  Table  1. 
Roundoff  noise  gain  factors  of  the  direct-form  cascade 
filters  are  also  given.  As  can  be  seen  from  the  table, 
better  roundoff  error  characteristics  are  obtained  by  using 
constrained  filter  design.  Roundoff  noise  gain  factors  of 
basic  and  minimum-noise  lattice  filters  are  given  in 
Table  2.  Note  that  the  roundoff  errors  of  minimum-noise 
lattice  filters  decrease  as  the  pipelining  stage  increases. 


Table  1  Roundoff  noise  gains  of  Normalized  laaice  and 
Direct-form  cascade  filters 


Normalized  lattice 

Cascade 

Constrained 

Look-ahead 

Look-ahead 

nonpipe. 

5.5246 

5.8544 

19874.1 

2-staee 

8.9003 

8.1268 

20211.9 

4-staee 

8.3172 

16.0518 

20227.5 

S-stage 

16.0489 

32.0286 

20340.6 

Table  2  Roimdoff  noise  gains  of  Noimalized 
and  Minimum-noise  lattice  filters _ 


Basic  lattice 

Minimum-noise 

Constrained 

Constrained 

nonpipe. 

4.5225 

2.1524 

2-staze 

9.0016 

2.0216 

4-suee 

8.5022 

1.1348 

8-suee 

16.1041 

0.2815 

The  number  of  multipliers  and  adders  needed  for  the 
implementation  of  these  filters  by  normalized  lattice 
filters  are  compared  in  Table  3.  Note  that  the  4th  order 
niter  by  constrained  design  could  be  implemented  with 
one  complex  pole  pair,  and  requires  fewer  multipliers 
than  all  other  filters.  The  table  shows  that  the  lattice 
filters  obtained  by  using  the  constrained  filter  design 
method  require  less  hardware  complexity  than  those 
obtained  using  look-ahead  on  the  nonpipelined  filter. 
Table  4  shows  the  number  of  multipliers  and  adders 
needed  for  the  implementation  of  basic  lattice  and 
minimum-noise  lattice  filters  obtained  by  using  the 
constrained  filter  design  method. 


Table  3  Number  of  Multipliers  and  adders  of  Normalized 
lattice  and  Direct-form  cascade  filters _ 


Normalized  latice 

Cascade 

Constrained 

Look-ahead 

Look-ahead 

Multi 

Adda 

Multi 

Adda 

Multi 

Adda 

nonpipe. 

19 

11 

19 

11 

5 

8 

2-5tage 

23 

15 

23 

15 

9 

12 

4-stage 

15 

11 

31 

23 

13 

16 

8-s(age 

23 

19 

47 

39 

17 

20 

Table  4  Number  of  Multipliers  and  Adders 
of  Basic  and  Minimum-noise  lattice  filters 


Basic  lattice 
(Constrained) 

Minimum-noise 

(Constrained) 

Multi. 

Adda 

Multi. 

Adda 

nonpipe. 

13 

12 

22 

11 

2-stage 

17 

16 

30 

15 

4-stage 

13 

12 

22 

11 

8-stage 

21 

20 

38 

19 

6.  Conclusions 

In  this  paper,  it  was  shown  that  the  lattice  filter 
structures  in  [1],  [2]  and  [4]  are  suitable  for  pipelining. 
By  using  constraint  filter  design  and  the  lattice  filter 
design  procedures,  we  designed  pipelined  high-speed 
lattice  filters  with  extremely  small  roundoff  errors.  The 
hardware  complexity  of  these  structures  is  also  much  less 
than  the  pipelined  lattice  filters  obtained  by  using 
scaaered  look-ahead  on  the  nonpipelined  filters. 

As  for  further  research  work,  the  development  of  new 
pipelined  lattice  filter  structures  is  also  of  interest.  Wave 
digital  filters  also  possess  very  good  finite  word  length 


properties.  Design  of  pipelined  wave  digital  filters  needs 
to  be  addressed  in  future.  In  [13]  -  [14],  most  significant 
bit  first  redundant  arithmetic  has  been  used  to  pipeline 
recursive  and  wave  digital  filters.  Our  pipelined  filters  can 
be  used  in  combination  with  the  approaches  presented  in 
[13]  -  [14]  also. 
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